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ABSTRACT 

We report the discovery of a system of unusual Hi filaments which appear 
to be associated with molecular clouds in the Perseus spiral arm of our Galaxy. 
We investigate the hypothesis that this system is the result of a directed flow of 
dissociated gas from clouds trapped within an extended wind flow from massive 
stars. The Wolf-Rayet star WR 5 and the Olb(f) star HD 17603 are identified 
as candidate driving sources. However, an examination of this hypothesis within 
the context of the theory of mass- loaded winds shows that these two stars alone 
can not account for the energetics and kinematics of the required spherically 
symmetric wind flow. Unless the apparent association between Hi, molecular 
gas, and stars is an accidental one, we suggest that other as-yet unidentified 
stars must have contributed to driving the filaments. 

Subject headings: ISM: bubbles — ISM: clouds — ISM: kinematics and dynamics 
— stars: individual (WR5, HD 17603) — stars: winds, outflows 
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1. Introduction 

Massive O-type stars, including the progenitors of Wolf-Rayet stars, inject enormous 
quantities of energy into the interstellar medium (ISM) in the form of dissociating and 
ionizing radiation and winds. The radiation fields create circumstellar /interstellar Hi and 
Hn regions into which winds blow, creating large-scale swept-up expanding shells of gas 
known as stellar wind bubbles (Dyson and de Vries 1972; Koo and McKee 1992a,b). There 
are numerous analytical and numerical studies of the formation, structure, and evolution 
of stellar wind bubbles (SWBs), but there are few convincing observations of Hi bubbles 
around Wolf-Rayet stars and massive O-stars. The Vllth catalogue of galactic Wolf-Rayet 
stars indicates that only « 15% are associated with Hi bubbles (van der Hucht 2001). 
Surveys are incomplete, but many of these proposed bubbles are not securely established. 

The interstellar environment of massive stars is observed to be extremely complex and 
is likely to be the combined result of inhomogeneous initial conditions and the influence of 
the stars. Our understanding of the physics of astrophysical flows in inhomogeneous media 
is incomplete (Hartquist and Dyson 1993). In particular the effects of conduction, diffusion, 
turbulence, magnetism, and of mixing and mass-loading processes are all poorly understood. 
Therefore, in the observational study of SWBs and the environment of massive stars it is 
important to remain open to the possibility of phenomena not predicted by standard SWB 
models. 

In this paper we present one possible example of such a phenomenon. We investigate 
the nature of a large-scale system of neutral atomic hydrogen filaments which apparently 
emanate from compact molecular clouds located in or near the Perseus arm of our Galaxy. 
We investigate the possibility that these Hi "tails" are the result of molecular cloud mass- 
loading within the wind flow of massive stars. 

The observations and data processing are briefly outlined in the next section. Section 3 
describes the observational results. In section 4 we consider possible driving sources, and in 
section 5 we take a closer look at the environment of the two main candidates. The discussion 
follows (§6), in which we first consider the kinematics and timescales involved (§6.1), and 
then present an analysis within the framework of a mass-loaded flow model (§6.2). Our 
conclusions follow in section 7. 



2. Observations and Data Processing 

The data discussed in this paper were obtained as part of the Canadian Galactic Plane 
Survey, a long-term multi-wavelength imaging survey of the ISM of the northern Galactic 
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plane (Taylor et al. 2003). Table 1 summarizes the main parameters of the relevant CGPS 
data. The principal observational components of the CGPS are arcminute resolution mosaics 
of 21 cm (in all Stokes parameters) and 74 cm (right hand circular polarization) radio 
continuum, and 21 cm Hi line emission, all observed with the Synthesis Telescope (Landecker 
et al. 2000) at the Dominion Radio Astrophysical Observatory of the Herzberg Institute of 
Astrophysics. 

Complementary data produced as part of the CGPS include mosaics of dust emission 
created using HIRES-processed IRAS data (Cao et al. 1997; Kerton and Martin 2000) and 
of 12 CO(J =1-0) emission from the FCRAO Outer Galaxy Survey (Heyer et al. 1998). 
The CGPS CO mosaics are the result of a reprocessing of the original OGS data, wherein 
correlated noise between pixels caused by the observing method was suppressed, and faint 
contaminating emission found in the reference positions of the survey was removed (Brunt 
et al. 2000). The reprocessed data were resampled and projected into the CGPS mosaic 
system. 



3. Results 

3.1. Neutral Atomic Hydrogen 

A mean brightness temperature Hi image, averaged between LSR radial velocities —53.4 
and —66.6 km s _1 , is shown in Figure 1. Prominent in this image is a linear Hi structure 
1?8 long, running at an angle ~ 20° north of west relative to the Galactic plane. In the 
extreme southeast this structure is a relatively narrow (~ 2' wide) filament beginning near 
(I, b) = (138?7, -2?0). Near (138?2, -1?9), the narrow filament is joined by a broader (~ 0?5) 
structure which runs parallel and to the south of it. This structure continues to broaden 
towards the northwest, reaching a maximum width of ~ 0?8. Its northwest extremity lies 
near (136?7, — 1?0), where it blends into larger scale Hi structures in the Galactic plane. 
The approximate centre of the linear structure is at (/, b) ps (137?7, — 1?6); we designate it 
HI G137.7— 1.6— 60, and for brevity refer to it as HI G137 in this paper. 

The radial velocity width of HI G137 is 15 km s -1 , and there is a distinct velocity 
gradient along its length: velocities become more negative towards the northwest. This 
velocity gradient clearly separates HI G137 from the Galactic Hi emission to the northwest 
(Figure 2). Individual channel maps (Figure 3) hint at the existence of linear substructures, 
resulting in the overall impression that HI G137 is composed of several parallel filaments. 
We can detect no clear velocity structure in the direction transverse to the length of the 
filaments: this may be due to a combination of sensitivity and angular resolution limitations 
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for individual filaments and confusion between them. 

We are unable to make a direct measurement of the distance to HI G137. There are 
no known Hn regions in the direction of the filaments, and there are large discrepancies 
between kinematic and photometric distances in the Galactic plane in this region in the 
sense that kinematic distances can be greatly over-estimated compared to the latter (Foster 
and Routledge 2003). 

3.2. Carbon Monoxide 

The CGPS CO data were searched for molecular clouds having positional and/or velocity 
coincidences with HI G137. We found no CO emission near HI G137 within the velocity range 
of the Hi structure. Widening the velocity range of our search revealed an irregular chain 
of CO clouds in the southeast of HI G137 at -44 < v < -50 km s" 1 extending ~ 1?4 
in a direction perpendicular to the length of the Hi structure. Of these, the cloud with 
the brightest peak CO line brightness (T R * = 4.0 K), which we dub CO G138. 1-1.9-48 
(CO G138 for brevity), lies at the convergence point of the broad-structured component of 
HI G137 (Figure 4). The velocity difference between CO G138 and HI G137 ranges between 
5 and 20 km s -1 . The northern narrow filament of HI G137 appears to emerge from between 
a gap in the chain of CO clouds. Within this gap there is a faint (1.3 K) compact cloudlet 
at —47 km s _1 embedded within very faint (~ 0.5 K) diffuse CO emission. 

4. Search for Driving Sources 

The positional association between molecular clouds and the convergent point of the H I 
filaments suggests that HI G137 might be the result of a directed flow of dissociated gas from 
the molecular clouds. If HI G137 is associated with the CO clouds at ~ —48 km s _1 , its 
distance is D ~ 2 kpc, the approximate distance to the Perseus Arm at this longitude. The 
source(s) of the radiation and/or winds driving HI G137 will be at the same distance. The 
former would have to be located on the sky such that the H I filaments point towards them 
and the CO clumps lie between them and the filaments in a more or less collinear fashion. 
The SIMBAD database was used to search for O, B, and Wolf-Rayet type stars within 0?5 
of (l,b) = (138?5,-2?0) (Table 2). Of the resulting ten stars, the Wolf-Rayet star WR 5 
and the O-type supergiant HD 17603 currently have substantial stellar winds and ionizing 
radiation. These two therefore merit closest attention. 
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4.1. Wolf-Rayet 5 and HD 17603 

WR 5 (HD 17638) is a Wolf-Rayet star of spectral type WC6 (van der Hucht 2001). It 
has a reddening of Eb-v — 0.96 (Ay ~ 3) and a photometrically derived distance of 1.91 kpc 
(van der Hucht 2001). 

HD 17603 is an 07.5Ib(f) star (Walborn 1973) located 7'.5 northwest of WR 5. Stars of 
type 0(f) are thought to be precursors of Wolf-Rayet stars. The available UBV photometry 
for this star (Haug 1970) is consistent with a reddening E B _v = 0.93 (A v ~ 3). The 
absolute magnitudes of Of stars are not well determined: however assuming that this star 
has M v ~ —6.3 (Garmany et al. 1982), it lies at a distance D ~ 2.2 kpc. 

Both WR 5 and HD 17603 satisfy the positional criteria to be considered candidate 
driving sources. In addition, both are hot, luminous objects with strong winds. For both 
stars, the apparent alignment of Hi— CO— star is very good, and these stars are likely to 
be the most luminous and energetic of those we have identified. Finally, their approximate 
distances are consistent with that proposed for HI G137. Stellar wind parameters for both 
stars are listed in Table 3. 

A final consistency check on the stellar and CO cloud distances can be made by inte- 
grating the Hi line profiles towards the stars from zero velocity to the velocity where an 
Hi column corresponding to E B -v ~ 1 is found. For the WR 5 and HD 17603 sight lines, 
reddening implies an Hi column of N ~ 4.5 x 10 21 cm -2 . Our Hi data shows that this 
column density is reached when the line profiles are integrated from zero velocity out to 
~ —50 km s _1 , consistent with the radial velocity of CO G138. 

4.2. The B Stars 

In addition to the Wolf-Rayet and Olb(f) star discussed above, there are eight B stars 
within ~ 0?5 of the tip of HI G137. Among them is EO Per, a rapidly varying B0 supergiant. 
A B0 supergiant has a mass of roughly 30 M and so was approximately type 07.5 while 
on the main sequence. According to Howarth and Prinja (1989), such a star has a wind 
terminal velocity of 2300 km s -1 , slightly lower than the lower bound quoted for both WR 5 
and HD 17603, and a mass loss rate of 10 -6 ' 7 M yr" 1 , an order of magnitude smaller than 
the lower bound for the other two stars. At present this star is likely to have a higher mass 
loss rate but much lower wind terminal velocity than WR 5 and HD 17603, as well as a lower 
ionizing flux. 

The remaining seven candidates have very limited observational data upon which to 
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base determination of their type and distance. If we assume they are located at D ss 2 kpc 
and lie behind extinction of Ay ~ 3, then all these stars have photometry roughly consistent 
with early B type stars on the main sequence, although we cannot rule out the possibility 
that some are actually O type stars. 

5. The Interstellar Environment of WR 5 and HD 17603 

Based on low resolution (9') observations, (Arnal 1992) proposed that WR 5 is associated 
with an ovoidal Hi deficiency, interpreted as a cavity, at velocities around —14 km s _1 . 
We note that integrating out to —14 km s _1 in the Hi line yields a column density only 
~ 1.5 x 10 21 cm' 2 , corresponding to Eb-v ~ 0.3, which is too small to be consistent with 
the observed reddening to the two stars. WR 5 and HD 17603 are well beyond the distance 
corresponding to that gas velocity and thus cannot be associated with the ovoidal structure. 
Our high resolution (1') observations show that the Hi distribution is highly structured in 
both morphology and velocity, and the "cavity" is by no means clearly distinguishable. In 
addition, regions of depressed H I brightness within this "cavity" correspond in position and 
velocity with much of the CO emission seen around —14 km s _1 , which suggests that the 
relative lack of H I emission is due at least in part to the hydrogen being in molecular rather 
than atomic form. 

For the velocity range within which we associate WR 5 and HD 17603 (—53.5 < v < 
—42.0 km s _1 ), our Hi data do not show any convincing evidence for an Hi shell surrounding 
a cavity. However, it is possible that a gas shell does exist but is not clearly visible in the 
H I observations because of its low contrast or because parts of it are molecular. The clear 
anti-correlation between the CO clouds and the Hi brightness distribution in the —53.5 to 
—42.0 km s _1 velocity range and the relative lack of CO emission within ~ 40' of the stars 
lend some support to the latter possibility. 

Marston (1996) cited evidence from IRAS Skyflux images for a l?35-diameter dust shell 
around WR 5. Such a shell is not obvious in the HIRES processed 12, 25, 60, and 100 /zm 
mosaics. To investigate further, we used the 60 and 100 /im data to produce dust temperature 
and optical depth images. The scattered dust concentrations interpreted by Marston as a 
shell correlate much better with the velocity-integrated CO data than with the velocity- 
integrated Hi data, some with CO emission in the velocity range around —14 km s _1 and 
others with the CO around —48 km s -1 . This suggests that what was interpreted as a shell 
is in fact the result of projection effects. 

We inspected Ha data for this region from the Virginia Tech Spectral-line Survey (Den- 
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nison et al. 1998) and the Wisconsin H-Alpha Mapper (Haffner et al. 2003) Northern Sky- 
Survey. The VTSS has an angular resolution of ~ 1(6 but does not provide velocity infor- 
mation. The WHAM Ha data has a velocity resolution of 12 km s _1 but a lower angular 
resolution ~ 1°. The VTSS data suggest that WR 5 and HD 17603 both lie within a 
region relatively devoid of Ha emission. The morphology of this "void" and thus its pos- 
sible association with the stars is uncertain because a comparison of the VTSS image with 
velocity-integrated CO emission suggests that absorption by the dust in molecular clouds 
strongly influences the observed Ha brightness distribution. 

We conclude that there is at present no conclusive evidence for an H I or dust shell/cavity 
structure around WR 5 or HD 17603, although the stars appear to lie in regions relatively free 
of molecular, neutral atomic, and ionized gas. Any gas shell surrounding this region would 
have to be of low contrast or be rather inhomogeneous in its composition and structure: the 
latter implies a "leaky" shell within which a hot wind bubble might not be maintained. 

6. Discussion 

The morphological relationships between the Hi tails, the CO clouds, and the stars 
raises the question of whether all these are physically related. The possibility of a simple 
coincidence cannot be ruled out, but an explanation for the morphology and velocity struc- 
ture of the H I tails would still be needed. In the remainder of this paper we investigate the 
possibility that winds from massive stars could be responsible for driving long narrow H I 
tails from compact molecular clouds in their vicinity. 

A viable physical model of HI G137 must be consistent with all of the observational data 
we have presented above. Obviously, we need a stellar mechanism for converting molecular 
to neutral atomic hydrogen and accelerating the Hi away from the molecular clumps in 
long narrow tails at a velocity of ~ 5 — 20 km s -1 . WR 5 and/or HD 17603 must be able to 
supply momentum and energy at the rates required over a time comparable to the dynamical 
timescale. The molecular clumps must also be able to survive within the environment of the 
star(s) for this time. 

One question that immediately arises is the association of Hi tails with only some of 
the CO clouds in the —44 < v < —50 km s" 1 velocity range. Standard Galactic rotation 
curves predict a slope of the distance— velocity relationship at this location in the Galaxy 
of ~ 150 pc (km s -1 ) -1 . The above velocity range thus potentially covers a line of sight 
range of ~ 1 kpc and the velocity resolution of our CO data (0.98 km s~ x ) corresponds to 
~ 150 pc. The empirical rotation curve of Brand and Blitz (1993) suggests the slope may be 
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even steeper. Given the potentially very wide range of distances for the CO clouds, only a 
small subset of the observed clouds may be within the range of influence of the tail-driving 
source(s). Additionally, Hi flows from other CO clouds may be subject to confusion if they 
have relatively low radial velocities. 

6.1. Kinematics and Timescales 

Determining the mass of extended Hi structures is always subject to wide uncertain- 
ties, not least because of confusion with unrelated emission. This is particularly acute for 
structures having large velocity dispersions. We determined the mass of HI G137 by first 
removing an estimate of the background Hi emission on a channel-by-channel basis. This 
background was calculated by a linear interpolation anchored to velocity channels bracketing 
the velocity range of HI G137 which were judged to be free of emission from the latter. The 
background subtracted cube was then integrated over a polygonal region enclosing HI G137. 
The solid angle of this region was 1.1 square degrees, and the mean Hi column density 
derived was N — 1.2 x 10 20 cm~ 2 . The total Hi mass is 330D 2 M and the total mass, 
assuming a hydrogen mass fraction of 0.77, is M = 425D 2 M . The uncertainty in the mass 
estimate is at least 25%. The physical length of HI G137 projected onto the plane of the sky 
is ~ 30D pc. 

Assuming D = 2 kpc for HI G137, this yields a total mass of the structure of M — 
1700 M = 3.38 x 10 36 g. Assuming a systemic radial velocity of the mass reservoir of 
—47 km s _1 , the linear momentum contained in HI G137 is P = Pq/ cos(i) = 1.05 x 
10 42 .D 2 / cos(i) = 4.20 x 10 42 / cos(z) g cm s _1 , where we have assumed that the space ve- 
locity of the flow can be deprojected from the radial velocity by using an inclination angle to 
the line of sight %. The mass-weighted mean outflow velocity is then 12.4/ cos(i) km s" 1 , and 
the kinetic energy is E = P /cos 2 (i) = 6.67 x 10 47 (P/ cos(i)) 2 = 2.67 x 10 48 /cos 2 (?) erg. 
The driving source (s) presumably have injected momentum and energy isotropically into the 
ISM. The observed opening angle of HI G137 is 9 20°: assuming that the three dimensional 
structure is conical and viewed at an inclination i, the true opening angle is 0sin(i) , and 
HI G137 subtends a fraction / = (1— cos(0sin(i)/2))/2 of a sphere. Therefore, / < 7.6xl0~ 3 
(< 0.1 ster), the maximum value occurring when % = 90° (when the outflow axis lies in the 
plane of the sky). Over the lifetime of HI G137 (discussed below), the driving source(s) have 
supplied a total (isotropic) momentum and kinetic energy of, respectively, P* = Pj {ftp) and 
P* = E/(feE), where the e factors represent the efficiency of the transfer of momentum and 
energy to the Hi. The minimum value of P* and P* is obtained by minimizing the following 
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two expressions with respect to i: 

e P P, 1 

P / cos(i) 

e E E* = 1 
P /cos 2 (i) 

For 9 = 20°, the minimum values are P* imi „ = 342P /e P = 1.43 x 10 45 /e P g cm s" 1 
(when i ps 55°), and E* tmin = 526E /e E = 1.40 x 10 51 /e E erg (when i m 45°). A simple 
conical outflow model yields for the flow timescale r = tq/ tan(i) = 2.36 x 10 6 P/ tan(i) 
4.72 x 10 6 / tan(i) yr, where we have used the mass- weighted mean outflow radial velocity 
12.4 km s _1 and the mass-weighted mean angular extent 1?72. We have noted above the 
presence of a strong uniform radial velocity gradient along the H I tails. If this were a mani- 
festation of velocity-sorting rather than continuous acceleration, r might be more accurately 
estimated by using the maximum outflow radial velocity observed (~ 20 km s _1 ) rather 
than the mean velocity, resulting in a factor of 2 decrease in r. We will argue below that 
the velocity gradient is more likely to be the result of a steady acceleration of the tails, in 
which case the use of the mean velocity is approximately correct. 

Using the dynamical timescale we can estimate the rate at which the driving source(s) 
deliver (isotropically) momentum and kinetic energy, respectively P* and E*. The minimum 
values of P* and P* are obtained by minimizing the two following expressions with respect 
to v. 

epP*r sin(i) 2sin(i) 




e£p*To sin(i) 2sin(i) 

P /cos 3 (*) ~ (l- cos (^)) cos 3 (*) 

We find P*, m j n = 341P /(epT ) = 9.6 x 10 30 /ep g cm s~ 2 (when i = 35°), and E^, )min = 
405Po/(e£;To) = 7.2 x 10 36 /es erg s _1 (when i = 30°). Figure 5 shows plots of the dependence 
of P*, P*, E*, and P* on i for the simple conical flow model. The uncertainty in the 
values of the kinematic and energetic parameters is dominated by our lack of knowledge of 
the inclination angle, affecting in particular the value of the kinematic timescale. In our 
discussions below we will assume that % most likely lies within the range 30° < i < 60°. Our 
estimates of the physical parameters of the H I tails are summarized in Table 4. 
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6.2. Mass-loaded Flow Model 

The lowest velocity shift discernable between the CO and Hi, limited by background 
confusion, is ~ 5 km s _1 , which is roughly comparable to the sound speed in neutral atomic 
hydrogen. The systematic radial velocity gradient of HI G137 could be interpreted in terms 
of curvature along the line of sight or of expansion. However, if HI G137 is the result of a flow 
from the CO clouds, it is reasonable to assume that the Hi filaments are radially directed 
away from the driving source and are thus approximately linear. The velocity gradient would 
then reflect a continuous acceleration of gas along the filaments, leading to the picture of a 
flow embedded within an extended wind region. This picture of long flowing tails from dense 
clouds driven by an external wind is similar to the situation described analytically by Dyson 
et al. (1993) and numerically by Falle et al. (2002) of the quasi-steady hydrodynamic mixing 
and entrainment of dense clump gas by a diffuse subsonic or transsonic wind. We propose 
that the neutral atomic tails are created in the dissociation and entrainment of molecular 
clouds enveloped within a strong stellar wind or winds. The wind ablates and entrains cloud 
gas and accelerates it away from the cloud in long narrow tails as a mass-loaded flow. Given 
the positional associations noted above, the source of the wind(s) is likely to be in the vicinity 
of WR 5 and HD 17603. 

We note that it is very unlikely that these Hi filaments are cometary tails caused by 
shadowing of radiation in a low-density Hll region. The flaring of the tails away from the 
axis of symmetry with increasing distance from the apex is inconsistent with H I observations 
(Moriarty-Schieven et al. 1996; Heyer et al. 1996) and models (Bertoldi and McKee 1990; 
Canto et al. 1998; Pavlakis et al. 2001) of cometary clouds. The dynamical timescale of 
HI G137, a few million years, is far longer than the timescale for complete photoionization of 
the tails (Pavlakis et al. 2001). This leaves the wind-driven model, but in order to establish 
the plausibility of this model there are several questions which must be addressed. Can the 
molecular clouds survive a few million years inside a stellar wind region? What must be 
the conditions inside this region in order to account for the observed cloud mass loss rate 
and entrained gas acceleration, and how realistic are they? Are the energetics plausible? 
What are the mechanisms by which the entrained gas may be dissociated? We address these 
questions in the following sections. 



6.2.1. Cloud Survival 

For the proposed model to be feasible, the molecular clouds must survive their initial 
impact with the stellar wind, before settling down to a quasi-steady entrainment phase 
within the wind. Klein et al. (1994) have performed an analytic and numerical study of 
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the interaction between a fast, steady, planar shock and a dense isothermal, nonmagnetic, 
spherical, cold cloud. More recently, Poludnenko et al. (2002) have extended numerical 
studies to the case of shocks encountering multiple clouds in which cloud-cloud interactions 
can oocur. As summarized in both these papers (we follow the formulation of Poludnenko 
et al. (2002) in this paper), a key timescale is the "cloud-crushing time" t cc , which is the 
time required for the internally transmitted shock to cross the cloud, t cc = 2ao/v cs . Here ao 
is the initial cloud radius and v cs , the forward velocity of the shock through the cloud, is 



where _F cl is the ratio of the pressure just behind the cloud shock to that at the stagnation 
point and F st is the ratio of the pressure at the stagnation point to that in the far upstream 
wind. Both these quantities have values of roughly unity (Klein et al. 1994; Poludnenko 
et al. 2002). The velocity of the external shock is v s and the density contrast x — Pco/pw 
is the ratio of the initial cloud density to that of the external medium. In the numerical 
simulations of Klein et al. (1994) and Poludnenko et al. (2002), the cloud destruction time 
tdest (defined as the time at which the mass of the surviving cloud core is reduced to ~ 0.5 
of its initial value) was typically found to be tdest ~ 2t cc . The numerical simulations of Mac 
Low et al. (1994) suggest values of tdest ~ 4t cc with the inclusion of dynamically significant 
magnetic fields. 

The molecular clouds at the apex of the Hi tails have a total mass ~ 440 M . Their 
current radii are a ~ 2 pc and they have H2 number densities % 2 ~ 150 cm -3 . At their 
present stage of evolution, the clouds have gravitational virial parameters a = 5a 2 a/Gm c ~ 
80-150, indicating that they are not bound by their self-gravity (<r is the observed cloud 
velocity dispersion, ~ 3.3 km s" 1 , and m c is the cloud mass). If we assume that the total 
mass of the H I tails was originally in the molecular clouds, we find that the clouds have been 
reduced from their initial mass by a factor of about 5: according to the formal definition of 
t dest above, these clouds would be classified as largely destroyed. 

If we set tdest = C^cc, the cloud destruction time is 



For shock velocities v s ~ 100 km s" 1 and an initial density contrast x > 100 (with a Q a few 
parsec and 2 < ( < 4), the clouds will survive for times of order a few 10 6 yr, long enough for 
them to become enveloped by the wind. We conclude that the molecular clouds can survive 
the passage of the initial shock, after which they are gradually destroyed while mass-loading 
the surrounding flow. 
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6.2.2. Conditions in the Mass-loading Region 

The mass-loading of flows by embedded clouds has been discussed in a series of papers 
by Dyson and collaborators (Hartquist et al. 1986; Dyson and Hartquist 1987; Charnley et 
al. 1988; Arthur et al. 1993; Dyson et al. 1993; Arthur et al. 1994; Dyson and Hartquist 
1994). For an initially pressure-confined cloud which becomes immersed in a wind flow, the 
pressure gradients across the cloud surface cause a lateral expansion of the cloud which mixes 
cloud gas into the wind. Hartquist et al. (1986) show that the mass loss rate from the cloud 
is 

m c = aM^(m c C c ) 2 ^( Pw v w f^ for M < 1 
m c = a(m c C c ) 2 / 3 (p w v w y/ 3 for M > 1 

where C c is the isothermal sound speed in the cloud, p w and v w are the wind density and 
velocity, and M = v w /C w is the Mach number in the wind. In this expression, a is a constant 
of order unity for which Hartquist et al. (1986) provisionally assign a value 1/2. Assuming 
that the wind and cloud parameters (except m c ) are roughly constant over the cloud lifetime, 
the cloud mass loss rate averaged over the entire cloud lifetime is (m c ) m m c0 /3, where m c0 
is the maximum mass loss rate corresponding to the initial mass m c0 of the cloud. Therefore, 

i™4 - 3 - 41 >< iqi8m4/3 ( — y /3 ( ^tV /3 (^y /s (^rV /3 

g s- 1 V M ©/ \ km s / Vcm" 3 / Vkms- 1 / 

where n w is the number density of the cloud (we have assumed a mean mass per particle of 
H = 2.17 x 1CT 24 g). 

For simplicity, we assume that the initial cloud mass (the current molecular cloud mass 
plus the mass in the H I tails) was contained within a single cloud, thus m c0 ~ 2100 M . The 
relative fraction of the initial mass in the tails and the clouds indicates that r ps 0.4^, where 
t H is the cloud lifetime predicted by integrating the cloud mass loss rate. We thus appear to 
be observing the clouds near the half-way point in time of their destruction process. The mass 
loss rate averaged over this time period (the dynamical timescale r) is (m c (r)) ~ 2m c o/3, 
roughly twice (m c ). Using our value of m c0 , we obtain 

(m c (r)) 21 / C c \ 2/3 / n w \V3 / Vw \V3 



g s 1 



\ km s 7 Vcm" - 3 ) (kms- 1 ) 



where we have assumed v w w C w (i.e. M 1). We have observed a value (m c (r)) 1.3 
- 3.9xl0 22 g s^ 1 (for 30° < % < 60°), which implies that we can match the theoretically 
predicted mass loss rate if 

1.5 x 10 3 < ( -^—) 2 < 4.2 x 10 4 . 
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In the theory of Hartquist et al. (1986), the isothermal sound speed C c is better identified 
as the characteristic speed at which the dense cloud can supply mass to the wind-cloud mixing 
region. Thus, we will assume that the effective value of C c is a rs 3.3 km s -1 . Then, the 
above condition becomes 



140 < — — — r < 3800 
~ cm" 3 km s- 1 ~ 

and has a value of ~ 730 when % = 45°. The wind velocity experienced by the dense 
clouds cannot be less than the maximum velocity observed in the Hi tails (~ 25-40 km s _1 , 
depending on i). The maximum possible wind speed is however much larger: the winds of 
massive main sequence stars can have terminal velocities ~ 3000 km s" 1 (cf. Table 3). If 
v w ~ 3000 km s -1 , the mass-loading model predicts 0.05 < n„,/cm~ 3 < 1.4. Such large 
values of n w rule out the possibility that the tails are produced in the mass-loading of clouds 
immersed in a freely-expanding bare stellar wind. Direct momentum-driving of the tails by 
a wind is also ruled out: the total isotropic linear momentum supplied by the wind of a 
40-85 M Q star during its main sequence lifetime of a few 10 6 yr is approximately 10 42 - 
10 43 g cm s" 1 (Schaerer et al. 1996), which is at least two orders of magnitude smaller than 
p 

1 *,min- 

The isotropic momentum flux predicted by the mass-loading model is 

P = 4nr 2 p w v w 2 = (9.2 - 83. x 10 24 )v w g cm s~ 2 , 

where substituting our adopted range of values of n w v w leaves only one factor of v w (the 
separation between the driving source(s) and the CO clouds is r = r /sin#, where we 
have assumed = 25 pc). Equating P with our observed value of 9.6 x 10 30 < P* < 
1.7 x 10 31 g cm s -2 , we see that wind speeds in the mass- loading region greater than v w ~ 
10 6 cm s -1 (~ 10 km s -1 ) would be sufficient to account for the momentum flux. Given that 
velocities in HI G137 may be as high as ~ 40 km s _1 , v w must be significantly higher than 
the limit of ~ 10 km s~ x we have obtained. We note that the predicted value of n w v w in the 
mass-loading theory is very sensitive to the value of the "mixing efficiency" a (n w v w oc a~ 3 ). 
The mixing efficiency may well be smaller (Canto and Raga 1991; Arthur and Lizano 1997). 
Lowering the value of the mixing efficiency by a factor of 2 or 3 would increase the required 
wind speeds in the mass-loading region by a factor of ~ 10, to ~ 100 km s -1 . Using this 
higher value, the mass-loading theory predicts a wind density in the mass-loading region of 
n w ~ 1-40 cm~ 3 . 

We conclude that the mass-loading scenario is workable in the context of the main 
sequence stellar winds from massive O stars if the fast low-density stellar wind is significantly 
loaded with mass and greatly slowed before interacting with the dense cloud. The low 
inferred wind velocity is consistent with the relative narrowness of the H I tails as predicted 
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theoretically (Dyson et al. 1993) and found numerically (Falle et al. 2002) for the case of mass 
loading by a mildly supersonic (transsonic) wind. A related point concerns the apparent 
complete mixing of mass injected into the wind before encountering the CO clouds: in 
contrast the H I filaments clearly do not completely mix over a large distance. Mixing may 
be greatly enhanced in a highly inhomogeneous medium having numerous clumps which are 
close enough togther for wind-driven interactions to occur (Poludnenko et al. 2002; Garcia- 
Arredondo et al. 2002). In contrast, the CO clouds associated with the Hi filaments are 
likely to be the dominant mass concentrations in their vicinity and so are not subject to 
mixing driven by cloud— cloud interactions. 

Our model still requires that the wind transfer very large amounts of kinetic energy 
to the Hi tails: several times 10 36 erg s -1 , which is the same order of magnitude as the 
mechanical luminosity of the main sequence wind of a massive O star. This is a severe 
constraint, since in the basic theory of adiabatic stellar wind bubbles, only a tiny fraction, 
~ 1%, of the energy is in kinetic form within the bubble, the bulk being in the thermal 
energy of hot shocked gas and in the kinetic energy of the swept-up shell. Therefore a very 
small fraction of the energy is available to drive mass-loading of clouds trapped inside. 

Much higher efficiencies of energy conversion into kinetic form are possible if the stellar 
wind is itself significantly mass- loaded. Arthur et al. (1994) have shown that mass-loading 
of an initially highly supersonic wind by mass sources distributed smoothly throughout 
the region around the star can reduce the Mach number of the wind without having the 
wind going through a global shock near the star. If the mass loading is sufficiently heavy 
(^load/^wind > 10, where mioad is the mass loaded into the wind and m win d is the mass 
of the wind injected by the star), a large fraction of the stellar wind energy (> 50%) will 
remain in the form of kinetic energy within the wind region (Pittard et al. 2001). There- 
fore, stellar winds which are heavily "pre-loaded" with mass may be able to provide the 
kinetic energy required to drive flows from larger-scale mass concentrations. It also may 
potentially explain the low wind velocities and high wind densities that we infer. In this 
regard, a scaled-up version of the scenario modelled numerically by Garcia- Arredondo et al. 
(2002), the mass- loading of the wind from the Trapezium cluster by propylds, may be rele- 
vant. Garcia- Arredondo et al. (2002) showed that mass-loading could have very significant 
dynamical effects on the wind, including increased wind density and low Mach numbers at 
large radii. However, even with a high efficiency of momentum transfer to the Hi filaments, 
the momentum input of the winds of WR 5 and HD 17603 integrated is too low by at least 
an order of magnitude. It would seem then that our scenario requires the wind input from 
a group of many more O stars in addition to these two. 
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We have not been able to identify such a grouping of stars. Since massive stars form 
in OB clusters, short-lived stars such as WR 5 and HD 17603 do not have time to travel 
far from their birthplaces (unless they are runaway stars ejected from a binary during the 
supernova explosion of the other component). Although it is not satisfactory to have to 
posit the existence of an unseen group of OB stars, the apparent close proximity of WR 5 
and HD 17603 to each other makes this idea more plausible. Given the effects of visual 
extinction (and possible anomalous extinction) in a direction towards a major spiral arm, 
and the difficulties of finding, classifying, and obtaining distance estimates, it is perhaps not 
unlikely that such a grouping of stars could have escaped detection. Given the dynamical 
timescale of the H I filaments, it is quite possible that the parent molecular cloud in which 
the group formed has been largely destroyed by now, and some of its members no longer 
exist after having gone supernova (we could speculate that supernova explosions might have 
contributed to the flow driving). Such a remnant grouping of stars might be very difficult 
to identify. 



6.2.3. Dissociation 

There are two plausible mechanisms for dissociating the molecular gas which forms 
the Hi tails: photodissociation and shock dissociation. Stars like WR 5 and HD 17603 
produce very large fluxes of ionizing photons (> 10 49 s" 1 ), and thus the ionized skins of 
the molecular clouds may be observable in Ha. The Ha brightness can be estimated using, 
e.g., equation 13 of Lopez-Martin et al. (2001), where we assume a cloud radius 2 pc, a 
ratio of ionized skin thickness to cloud radius of 0.1, and a source of ionizing photons ~ 
7.5 x 10 49 s _1 located 25 pc from the cloud. The latter was estimated using equation 1 
in Greiner et al. (1999). The resulting Ha surface brightness observed from the Earth, 
assuming no extinction, corresponds to a photon specific flux of ~ 3.6 x 10 7 s _1 cm -2 ster -1 . 
This is comparable to the surface brightnesses of the fainter photoevaporating cometary 
knots observed in the Helix nebula by O'Dell et al. (2000). If the clouds suffer a similar 
extinction as WR 5, the Ha emission will be further attenuated by a factor of ~ 10. The 
large (Ay ~ 3) and patchy visual obscuration may make detection of Ha emission from 
the clouds difficult (see discussion in section 5). There is no evidence for radio continuum 
emission from ionized hydrogen from the CGPS 1420 MHz continuum data. Any Hll region 
must thus be of very low radio surface brightness: a 100 pc path length through a 1 cm~ 3 
H II region would produce only a few tenths of a kelvin brightness temperature, which could 
not easily be distinguished from the general Galactic background of several kelvin. 

Hot massive stars produce photons capable of dissociating molecular hydrogen at rates 
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comparable to the ionizing photon rates (Diaz-Miller et al. 1998). A rough evaluation of 
the possibility that the tails are photodissociated can be made by comparing the photodis- 
sociating rate into the solid angle of the tails, pSpf, to the dissociation rate implied by 
our observations. Sp ~ 10 49 s _1 is a typical stellar photodissociating photon rate for a 
very massive O star, p 0.15 is the probability that an absorption will lead to dissociation 
(Diaz-Miller et al. 1998), and / is the fractional solid angle subtended by the tails. We find 
that the two rates are roughly equal at ~ 10 60 dissociations when integrated over a period 
of a few 10 6 yr. Photodissociation thus appears to be a plausible dissociation mechanism for 
the tails. 

The other possibility is that shocks driven into the molecular clouds by the stellar 
wind(s) may be capable of dissociating the H 2 . J-type shocks, in which the shock velocity 
v cs is sufficiently high to preclude the existence of a magnetic precursor, are the most effective 
at dissociating molecular gas. In a strong J-type shock with v cs > 50 km s" 1 moving into 
molecular gas with a pre-shock H 2 number density n ~ 100 cm -3 , about 90% of the H 2 will 
be dissociated (Hollenbach and McKee 1980). For \ ~ 10 2 an d a wind velocity v w of the 
order of a few hundred km s -1 , shock velocities will be of the order of 40-50 km s _1 , which 
is sufficient to result in very significant dissociation. The energy required to shock dissociate 
the molecular gas is roughly an order of magnitude smaller than that to accelerate it to 
the velocity of the tail. However, shock velocities of this magnitude may reduce the cloud 
destruction time to values < 10 6 yr, which is a problem for the scenario we are proposing. 

We conclude that photodissociation is the most workable mechanism for dissociating 
the molecular gas swept into the H I filaments. 

7. Summary and Conclusions 

We have discovered an unusual system of H I "tails" comprised of a number of linear or 
quasi-linear thin filaments which, if projected along their length, appear to have a common 
origin. The feature, which we dub HI G137.7— 1.6— 60, has dimensions of l?8x0?8x 15 km s -1 
(length x width x velocity), although its width varies along its length. 

Located at one end of the filaments is a chain of molecular clouds that lie at a distance 
of approximately 2 kpc. Despite the molecular material being displaced from the atomic by 
~ 10 — 30 km s -1 , we propose that the molecular clouds and Hi feature are located at the 
same distance and are dynamically related. The lack of a traditional head-tail morphology, 
and the lack of any evidence of ionization along the leading edge of the molecular clouds, 
lead us to reject the possibility that this system is simply a cometary cloud; instead, we 
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propose that the Hi feature is made up of material dissociated from the molecular clouds 
and entrained in the outflowing wind as a mass-loaded flow. 

Given the velocity difference between the atomic and molecular material, we estimate 
the dynamical age of the system to be on the order of a few 10 6 years. The velocity gradient 
observed along the H I feature is, in our model, explained by continuous acceleration of the 
entrained material by the stellar wind. 

A search for driving sources for the H I filaments identified the Wolf-Rayet star WR 5 and 
the Of star HD 17603 as candidates. Both these stars are located near the apparent origin 
of the filaments, are approximately at the same distance as the molecular clouds, and have 
main sequence lifetimes comparable to the dynamical age of the system. However, even with 
the high efficiency of momentum transfer predicted for a heavily mass-loaded wind, WR 5 
and HD 17603 are incapable on their own to drive the filaments: we are forced to speculate 
that other, unidentified stars and/or old supernova explosions must also have been involved. 
We cannot of course discount the possibility that the positional associations between Hi, 
CO, and stars are simply accidental. But, assuming that stellar winds (or explosions) have 
shaped the gaseous structures we have observed, we suggest that anisotropic Hi structures 
similar to HI G137 should be searched for in the environments around massive stars (and 
supernova remnants). 

We thank Tom Landecker for his support and encouragement, and the referee for com- 
ments that led to a number of key improvements in this paper. The Dominion Radio Astro- 
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with international partners, and is supported by grants from the Natural Sciences and Engi- 
neering Research Council of Canada. Guest User, Canadian Astronomy Data Centre, which 
is operated by the Herzberg Institute of Astrophysics, National Research Council of Canada. 
Data from the Canadian Galactic Plane Survey is publicly available through the facilities of 
the Canadian Astronomy Data Centre (http://cadc.hia.nrc.ca). 
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Table 1. Parameters of Relevant CGPS Data 



Data Product 


Parameter 


Value 


Hi 21 cm line 


Spectral resolution 


1.32 km s" 1 




Channel width 


0.82 km s- 1 




Angular resolution 


64" 


12 CO(J = 1-0) line 


Spectral resolution 


0.98 km s- 1 




Channel width 


0.82 km s" 1 




Angular resolution 


100" 
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Table 2. All O, B and Wolf-Rayet stars within 0.5° of the tip of HI G137. M 



Name 


Spectral type 


G.long. 


G.lat. 




Comments 


WR 5 


WC6 M 


138.87 


-2.15 






HD 17603 


07.5Ib(f) W 


138.77 


-2.08 






EO Per 


B01ab:e 


138.50 


-1.89 


1 


Listed RA/dec differ slightly 


CDS 314 


B 


138.50 


-1.89 


J 


for these two stars 


BD+56°722 


B 


138.57 


-2.14 






ALS 7516 


B 


138.80 


-2.08 






ALS 7523 


B 


138.87 


-2.04 






ZZ Per 


B 


138.84 


-2.19 






ALS 7453 


B 


138.22 


-2.32 






ALS 7470 


B 


138.16 


-1.75 







^ Compiled using the SIMBAD database. Unless otherwise noted the information 
in this table is that listed in SIMBAD. 

Mvan der Hucht (2001). 

WWalborn (1973). 
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Table 3. Wolf-Rayet 5 and HD 17603 



Parameter 


WR5 


HD 17603 


Spectral type 


WC6 [a] 


07.5Ib(f) M 


v oo 






M 


10" 4 - 3 M yr- 1 W 


10- 6 - 79 M yr- 1 W 


L w 


8.9 x 10 37 erg s" 1 




MS spectral type [0] 


05 


05 


MS Voo 


2800 km s" 1 


2800 km s- 1 


MS M 


lO" 5 - 7 M yr- 1 


10~ 5 - 7 M yr" 1 


MS L w 


4.5 x 10 36 erg s" 1 


4.5 x 10 36 erg s -1 


MS lifetime 


3.7 Myr 


3.7 Myr 


MS wind energy 


5 x 10 50 erg 


5 x 10 50 erg 



'°lThe main sequence spectral types were estimated from the 
evolutionary tracks by Maeder (1990). WR5 would have had 
a ZAMS mass of 40 M yr -1 or more which corresponds to a 
main sequence spectral type of 06.5 or earlier. The type 05 
was chosen to give a rough estimate of the energies involved. 
Masses quoted for 07.51 stars are on the order of 60 M yr -1 
implying a main sequence type of 05 or earlier. Again, the 
type 05 was chosen to give a rough estimate of the energies 
involved. Their stellar wind parameters were calculated from 
the empirical relations of Howarth and Prinja (1989). Their 
main sequence lifetime is from Chiosi et al. (1978). 

l^van der Hucht (2001). 

^Rochowicz and Niedzielski (1995). Values found in the lit- 
erature range from 1600 km s -1 (Koesterke and Hamann 1995) 
to 2800 km s" 1 (Torres et al. 1986). 

^Koesterke and Hamann (1995). 
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[^Walborn (1973). 

HBieging et al. (1989), based on comparison with other stars. 

[^Bieging et al. (1989), based on empirical relations by Gar- 
many & Conti (1984). 



-22 - 



Table 4. Physical parameters of the H I tails 



Parameter Notation Value 



Average Hi column density in HI G137 


N 


1.2 


X 


10 20 


cm 2 


Total mass of HI G137 


M 


1700 


M 




Flow time scale for HI G137 


T 


4.7 


X 


10 6 ■ 


tan(i) 


Linear momentum in HI G137 


P 


4.2 


X 


10 42 


cos(j) ° 


Kinetic energy in HI G137 


E 


2.7 


X 


10 48 


— ttt erg 


Minimum total kinetic energy from driving sources 


E 

- LJ *,mm 


1.4 


X 


10 51 


— erg 


Minimum total momentum from driving sources 


P 

1 *,min 


1.4 


X 


10 45 


— g cm s _1 


Minimum rate of kinetic energy delivery 


E 

- LJ *,min 


7.2 


X 


1Q 36 


— erg s" 1 


Minimum rate of momentum delivery 


p 

1 *,min 


9.6 


X 


1Q 30 


1 —2 

— g cm s 
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Fig. 1. — Mean Hi brightness temperature image of HI G137, averaged between LSR radial 
velocities —53.4 and —66.6 km s _1 . 

Fig. 2. — The velocity gradient along HI G137 is highlighted by the transition from red (Hi 
emission at —57.1 km s -1 ), to green (—60.4 km s" 1 ), to blue (—63.7 km s" 1 ) as one moves 
along along the tails. The width of the velocity bins is 3.3 km s _1 . 

Fig. 3. — Eight Hi channel maps of HI G137.7— 1.6— 60 showing filamentary substruc- 
tures. The positions of WR 5 and HD 17603 are indicated. Each velocity channels spans 
0.82 km s _1 . Greyscale runs from K (white) to 75 K (white). 

Fig. 4.— Contours of 12 CO(J = 1 - 0) at -46.81 km s" 1 at 0.75, 1.25, and 2.0 K superim- 
posed on a greyscale of the Hi at —58.35 km s _1 . The northern narrow Hi filament appears 
to emanate from a faint CO cloudlet within a gap in a chain of brighter CO clouds. The 
thicker southern part of the Hi tails appear to emanate from CO G138. 

Fig. 5. — Dependance of P* (— — ), E* (• • • •), P* (— — ), and (—• — •—) on the 

inclination angle i. The first two parameters are scaled by ^- and the last two by -^S 
where X is P or E. 
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